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1. Introduction 

The modern experiments provide a possibility to discover new physical effects caused by pres- 
ence of the strong (hadronic scale) magnetic field. At the Relativistic Heavy Ion Collider (RHIC) at 
the first moments (t ~ 1 fm/c) of noncentral collision the very strong (B ~ 10 15 T, \fe~B ~ 300 MeV) 
magnetic fields appear[[I],^|]. Such strong magnetic fields can be also created in ALICE experiment 
at LHC, at the Facility for Antiproton and Ion Research (FAIR) at GSI and in the experiment NICA 
in Dubna. The additional motivation for the study of the effects induced by the strong magnetic 
field could also come from the physics of the early Universe, where the strong fields (B ~ 10 16 T, 
\fe~B ~ 1 GeV) could have been produced after the electroweak phase transition[|3|]. Due to the 
nonperturbative nature of the effects we perform the calculations in the lattice gauge theory. We 
use quenched approximation and show that for some problems it provides rather reasonable values 
of the physical quantities. 

This work has been done analogously to the previous SU(2) lattice studies [|], ||, ^, ffl]. The list 
of considered effects induced by the magnetic field is the following. 

The strong magnetic field can enhance the chiral symmetry breaking. There are various models 
(see Sec.|3j) which predict the growing of the chiral condensate. 

The second effect is the chiral magnetization of the QCD vacuum. This effect has a paramag- 
netic nature. The vacuum magnetization is related to the nucleon magnetic moments[||] and other 
nonperturbative effects of hadrons^]. We calculate the magnetic susceptibility and other quantities 
in Sec.|. 

The quarks develop an electric dipole moment along the field due to the local fluctuations of 
the topological charge^. We study this effect in Sec.||. 

Finally, the fluctuations of the topological charge can be a source of the asymmetry between 
numbers of quarks with different chiralities created in heavy-ion collisions. The so called "event- 
by-event P- and CP- violation" [[I]] can be explained by this asymmetry and observed at RHIC. So, 
our aim is also to see any evidences of this effect in SU(3) lattice simulations, nevertheless they are 
similar to SU(2) lattice results^]]. 



2. Technical details 



We use the quenched SU (3) lattice gauge theory with tadpole-improved Luscher-Weisz action 
[10]. To generate the statistically independent gauge field configurations we use the Cabibbo- 
Marinari heat bath algorithm. The lattice size is 14 4 , and lattice spacing a = 0.105 fm. All 
observables we discuss later have a similar structure: (^Pi^P) for VEV of a single quantity or 
( x i'€?i x ¥ ^^2^) for dispersions or correlators. Here &, G\, &i are some operators in spinor and 
color space. These expectation values can be expressed through the sum over M low-lying 1 but 
non-zero eigenvalues iXt of the chirally invariant Dirac operator D (Neuberger's overlap Dirac 



'We believe that the IR quantities are insensitive to the UV cutoff realized by selecting some finite number of the 
eigenmodes||ll|] 
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operatorJT2|]): 



wt^Wk 



\k\<M 



iXk + m 



(2.1) 



and 



k, P 



(k^^jp^p) - {p\^\k){k\^ 2 \p) 
(iXk + m)(iX p + m) 



(2.2) 



where all spinor and color indices are contracted and we omit them for simplicity. The are 
defined by the equation 



Dy k = iX k xlf k , 



(2.3) 



where \j/k are the corresponding eigenf unctions and the uniform magnetic field F12 = #3 = B is 
introduced as described mM\. To perform calculations in the chiral limit one calculates the expres- 



sion (gjj) or ( |2.2[ ) for some non-zero m and averages it over all configurations of the gauge fields. 
Then one repeats the procedure for other quark masses m and extrapolates the VEV to m — > limit. 



3. Chiral condensate 



In this section we present our results for the chiral condensate 



£ = -(0|W|0), 



(3.1) 



as a function of the magnetic field B. The general tendency for £ to grow with B was already 



obtained in various models: in the chiral perturbation theory [ |13| , |14[ ] (£ oc B for weak fields, E °c 
B 3 / 2 for strong fields), in the Nambu-Jona-Lasinio model [15] (£oc B 1 ), in a confining deformation 
of the holographic Karch-Katz model jl6| ] (£ <=c B 2 ), in D3/D7 holographic system Jl7| ] (£ oc fi 3 / 2 
for low temperatures, £ <*= B for high temperatures) and in SU(2) lattice calculations [Q](r °< B). 
Here our aim is to see how the chiral condensate behaves in the SU (3) quenched gluodynamics. 



We use the Banks-Casher formula [18], which relates the condensate (3.1) with the density 
p (A ) of near-zero eigenvalues of the Dirac operator: 



lim 



V ' 



where V is the four- volume of the Euclidean space-time. The result is shown in Fig. |l(a) 
We perform the fit of the results by the following function: 



Z fit (B) 



, eB 

1 + 1 i* 



where £0 = £(0). The obtained fitting parameters are 



(3.2) 



(3.3) 



[(228±3)MeV]\ A B = (1.31 ±0.04) GeV, v = 1.57±0.23. 



(3.4) 
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Figure 1: Chiral condensate 



It is interesting to compare quantitatively the condensate obtained by the Banks-Casher for- 
mula and that one calculated by the expression (2.1) with & = 1. The result is shown in Fig. l(b)| . 
The value of the condensate in absence of the magnetic field equals £(0) = [(230±5)MeV] 
which is not so far away from the value, which can be estimated by the Gell-Mann-Oakes-Renner 



formula! 22]: 



1(0) 



F 2 m 2 
1 n' n n 



2(m u + md) 



[(240 ± 10) MeV} 2 



(3.5) 



4. Chiral magnetization and susceptibility 

In this section we calculate the quantity 



(4.1) 



1 



where a a p = — \y a ,Yp\ and %{F) is some coefficient of proportionality (susceptibility), which 
depends on the field strength. 

This quantity was introduced in[|8|] and can be used to estimate the spin polarization of the 
quarks in external magnetic field. The magnetization can be described by the dimensionless quan- 
tity jU = % ■ qB, so that 



OFff^) =/i(W>. 



(4.2) 



The expectation value (4T) can be calculated on the lattice by (2T) with & = a a p. The result 



is shown in Fig. 2(a) (here for comparison we also plot series for some finite quark mass). We can 
see, that the 12-component grows linearly with the field, which agrees with[^]. This allows us to 
find the chiral susceptibility ^(0) = Xo- After making a linear approximation (Wan^) = 0.f"eB, 
where 2 



1 



(4.3) 



in our simulation we calculate the magnetization of the d-quark condensate, thus q - 
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Figure 2: Expectation values of WdaffV and their square 



we obtain Q.^ = (172.3 ±0.5)MeV and 



AO 



-4.24 ±0.18 GeV " 



(4.4) 



This value fits well into the range of present theoretical estimations: the modern QCD sum rule 
calculations S (xf = -3. 15 ±0.3 GeV" 2 ) and [|24|](-2.85 ± 0.5 GeV" 2 ), the earlier ones ^ 



{-5.1GeV- z ) and [M (-4.4±0.4GeV" 2 ), in the instanton vacuum model [1271] (-4.32G<?V" 2 ), 



r-2\ 



from the analysis of the Dirac zero-mode in an instanton background [^] (—3.52 GeV 2 ), in 
dubbed quark-meson model [^] (-4.3 GeV" 2 ) and in the Nambu-Jona-Lasinio model [g] (-5.25 GeV~ 2 ). 
We also have to mention a well known analytic result obtained by OPE combined with the idea of 
pion dominance [30] and two holographic ones 32], but for us it seems difficult to compare 
these results with ours because of an ambiguity in the determination of the pion decay constant F n 
in the quenched approximation. 

Another interesting phenomenological quantity is the product of the chiral susceptibility % and 
the condensate ( X P X P)[^]. In our calculations it is equal to 



52 MeV, 



(4.5) 



while from the QCD sum rules one can estimate this quantity as approximately 50 MeV [23, 25 
26h, which is also close to our value. 



5. Electric dipole moment 

Another interesting effect due to the magnetic field is a quark local electric dipole moment 
along the field[]6|]. This quantity corresponds to the /0-components of the ([Q|): 



1,3 



(5.1) 



In the real CP-invariant vacuum the VEV of this quantity should be zero: (di(x)} = 0, that we 



actually see in our results (Fig.2(a^). At the same time the fluctuations of di(x) can be sufficiently 
strong. We measure VEV's (2.2) with 6\ = &2 = In the case of /0-components it corresponds 
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Figure 3: Fluctuations of the chirality and electromagnetic current/charge 



to dispersions of d. The result is shown in Fig. 2(b) , we see that the longitudinal fluctuations of the 
local dipole moment grow with the field strength, while transverse fluctuations are absent. Here 
and after we use the "IR" subscript to emphasize, that we subtract from the quantity its value at 
B = 0: 



{Y)m(B) = ^J d 4 x{Y(x)) B -^ J d 4 x{Y(x)) B=0 



(5.2) 



6. Some evidences of the Chiral Magnetic Effect 

The nontrivial topological structure of QCD is due to some nontrivial effects in the presence of 
the strong magnetic field. One example of a such effect is the Chiral Magnetic Effect (CME), which 
generates an electric current along the field in the presence of the nontrivial gluonic background^, 



21[ ]. This effect was probably been observed by the STAR collaboration at RHIC[]19[, [20| ] in heavy- 
ion collisions. A lattice evidence of the effect can be found in[f7|, [3^, 34]. Here we implement the 
same procedure for the SU (3) case and study the local chirality 



p 5 (x) = y(x)Y 5 V(x) = p L (x) - p R (x) 



and the electromagnetic current 



j^x)=^( x )r^(x). 



(6.1) 



(6.2) 



The expectation value of the first quantity can be computed by ( |2.1[ ) with 6 = y 5 and with G = 
for the second quantity. The both VEV's are zero, as expected, but the corresponding fluctuations 
obtained from (22) are finite and grow with the field strength (see Fig.||). 
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